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Fidelity decay is studied for quantum many-body systems with a dominant independent particle 
Hamiltonian resulting e.g. from a mean field theory with a weak two-body interaction. The diagonal 
terms of the interaction are included in the unperturbed Hamiltonian, while the off-diagonal terms 
constitute the perturbation that distorts the echo. We give the linear response solution for this 
problem in a random matrix framework. While the ensemble average shows no surprising behavior, 
we hnd that the typical ensemble member as represented by the median displays a very slow hdelity 
decay known as "freeze" . Numerical calculations confirm this result and show that the ground state 
even on average displays the freeze. This may contribute to explanation of the "unreasonable" 
success of mean field theories. 

PACS numbers: 03.65.Yz,05.30.Fk 
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I. INTRODUCTION 

Wigner [l[ proposed the use of random matrix the- 
ory, in particular of the classical ensembles [1, [1] to ex- 
plain statistical properties of nuclear spectra. Balian Q 
showed that this amounts to using a minimum infor- 
mation approach, where symmetry is the only charac- 
teristics taken into account. An essential element for 
the widespread success [5] resides in the ergodicity of 
the results obtained for classical ensembles. Yet it was 
soon noticed [1, 0] that the essential two-body character 
of the underlying interaction requires embedded ensem- 
bles, in particular the two-body random matrix ensemble 
(TBRE). While the original work carried the full weight 
of a three dimensional few nucleon system, soon abstract 
models with structureless fermions were introduced Q 
to understand more easily embedded ensembles in gen- 
eral and the TBRE in particular. Two problems beset 
these studies: one is the apparent non-ergodicity of these 
ensembles, the other is the difficulty to obtain any ana- 
lytical results (see Q for a review). After individual un- 
folding of each spectrum [1, 0| spectral statistics of the 
classical ensembles are recovered and this marginated the 
interest in embedded ensembles for some time. Interest 
was rekindled recently mainly by studies of interesting 
properties at the edge of the spectra concerning ground 
states [13, [m, ■ The question whether the ensemble 
is actually ergodic or not remains open but French 
has shown that it is somewhat academic, as even 
the spectral density has fluctuations of the order 1 / log N 
where TV is the dimension of the Hilbert space. This 
makes individual unfolding essential for any statistical 
analysis of spectra. The reason why individual unfolding 
does lead to the right answers is not understood to this 
day. The renewed interest in the TBRE has spread to 
various fields includirig applications @, [l4|, [l5| and ex- 
tensions to bosons d, [iBl ■ 



Fidelity analysis for the stability of quantum systems 
under perturbation, first proposed by Peres [l3|, has be- 
come very fashionable since the advent of quantum in- 
formation where it provides a standard criterion of 
stability; for a recent review see [l^. It seems reason- 
able to check how a small residual two-body interaction 
affects the stability of the solution of an independent par- 
ticle Hamiltonian. The diagonal part of this interaction 
does not affect the eigenstates and is usually included in 
the independent particle term, known as the mean field 
approximation. 

We focus on the validity of the mean field approxima- 
tion, given by the stability of the time-dependent solution 
of the mean field Hamiltonian which is measured by fi- 
delity. This is also the main physical motivation behind 
defining the random matrix model in the present paper. 

Fidelity has been calculated by the method of sta- 
tionary phase applied to time-depende nt p ropagator for 
semi-classical considerations [10, Hi], IH, [2J| and by super- 
symmetric techniques [2^ [25j for a random matrix model 
presented in ^261]. In both contexts most of the relevant 
results can be obtained by perturbative or linear response 
techniques, which actually do not require either model 
as a basis [19, 26, 27, 28, 29, 30]- Even in the regime 
where fidelity is not anymore close to one, a simple expo- 
nentiation of the lea ding term yields excellent results far 
beyond this regime [mill. Cr ossovers between semi- 
classical and perturbative considerations (the latter also 
known as Fermi Golden rule regime) have been discussed 
in [3l|, H^l ■ A drastic reduction of fidelity decay termed 
"fidelity freeze" occurs when the perturbation is strictly 
off diagonal (or more generally, when its time average 
is zero) (sol , [ssl [3^ . Such a situation can occur when 
perturbation breaks the antiunitary symmetry (like time- 
reversal) of the unperturbed Hamiltonian in an optimal 
way [30|, which can perhaps be experimentally realized 
by some magnetic interactions. 
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In case of high level fidelity freeze the dynamics is sta- 
ble for a very long time. The question is whether we can 
see this effect in the situation under consideration. On 
the one hand our system looks promising, due to the in- 
clusion of the diagonal part of the interaction in the "un- 
perturbed" independent particle Hamiltonian. On the 
other hand the occurrence of the freeze in random ma- 
trix models depends heavily on level repulsion implicit in 
a GOE or a GUE for the unperturbed system. Clearly 
the independent particle model does not have this repul- 
sion and indeed we can expect results similar to those 
of a random spectrum where absence of a diagonal per- 
turbation to lowest order only suppresses the quadratic 
decay 

Our analysis will show that the ensemble average does 
not display fidelity freeze. The typical ensemble member 
however does show this freeze. This is confirmed when 
we consider the median of fidelity rather than its average. 
This fact is important by itself as it shows that quantum 
evolution of the full problem will follow that of the in- 
dependent particle problem or mean field approximation 
for a very long time for most systems. It also is im- 
portant because it shows that RMT models can describe 
physical situations qualitatively, even when the average 
behavior is completely off. The last point can also be re- 
formulated as a necessity to look at the right quantities 
i.e. quantities whose distributions have no long tails. If 
we had replaced fidelity by distortion as introduced in an 
elastic problem [s^, then the average and median would 
coincide. This results because distortion is defined as 
the logarithm of one minus fidelity (36| . Indeed in other 
contexts taking averages of the logarithm is a common 
remedy to eliminate exaggerated effects of tails in distri- 
butions (37} . Furthermore we will see that for the ground 
state and the first excited state of the independent par- 
ticle model the freeze will occur even on average. 

After introducing basics about fidelity and the Hamil- 
tonian we will sec that under specific assumptions for the 
independent particle spectrum we can obtain ensemble 
averaged fidelity decay of a TBRE in the linear response 
regime. The result essentially yields the linear decay. Yet 
a numerical inspection of individual Hamiltonians in the 
ensemble shows the existence of fidelity freeze. A more 
careful numerical analysis shows that the median fidelity 
for the ensemble indeed displays the freeze. Thus we can 
consider it as typical. This result is quantitatively recon- 
firmed by showing that the ensemble averaged logarithm 
of one minus fidelity also displays the freeze. We relate 
this behavior directly to the existence of a gap in the 
nearest neighbor spacing distribution of the independent 
particle spectrum. Furthermore we find that the fidelity 
decay of the ground and first excited state of the inde- 
pendent particle Hamiltonian display the freeze even on 
average. Next we numerically check different options for 
the single particle spectrum to obtain a better feeling 
of the non-ergodicity, which plays a central role in this 
context. Finally we give some conclusions. 



II. FIDELITY AMPLITUDE 

The fidelity amplitude measures the overlap between 
two quantum states ipoit') and ip(t') evolving from the 
same initial state '0(0) but propagated with slightly dif- 
ferent Hamiltonians Hq and H = Hq + XV respectively 

fit') = {Mt'mt')) = {m\e"-''"''e-i''"\m)- (1) 

For small perturbation strength A it is convenient to 
use a linear response or Born expansion fiol . [27l . [28j . We 
shall follow the notation established in 26] as it is well 
adapted to RMT. In the interaction picture the wave 
function reads as x{f) = e's^"* Time t' can be 

replaced by dimensionless time t measured in units of 
the Heisenberg time tn = 2Trh/d where d denotes the 
average level spacing in the spectra of Hq which can and 
will be set to one. Therefore, the time will always be 
given in units of tn- Time evolution of a state x{t) up 
to the second order x^^-'(t) — X{t)x{0) in perturbation 
strength A can be expressed as 

X{t) = 1 - 27rAi [ dtiViiti) 
Jo 

- (27rA)2 f dh r dt2Vi{h)Vi{t2) (2) 
Jo Jo 

where Vi{t) is the abbreviation for the perturbation in 
the interaction picture Vi{t) = Q'i^'iHay^-2TTiHo _ piease 
note that the expansion in ([2|) is valid even for times 
much longer than the Heisenberg time [t ^ 1), if only 
the perturbation strength is small enough. 

The fidelity amplitude for an initial state |x(0)) = 
Xa\c») written as a superposition of independent par- 
ticle states, i.e. state eigenstates of the unperturbed 
Hamiltonian \a), then reads 

N 

fit) = (x(0)|X(t)|x(0)) = ^ X^0xlxp. (3) 

a, 13=0 

Later we consider ensemble averages. Then the linear 
term in ([2]) vanishes and the matrix element Xap is re- 
duced to contributions from the quadratic term only. 

III. THE HAMILTONIAN 

Although a simple RMT of full Gaussian matrices often 
yields impressively accurate results, it, strictly speaking, 
only applies to dynamical systems where all levels are 
coupled by the interaction. More realistic physical sys- 
tems such as nuclear or atomic shell models or quantum 
dots do not possess that property. The interaction in- 
volves two particles only. In the framework of RMT it 
should be described by Two-Body Random Ensembles 
(TBRE) 0, S m. We consider a system of M 

"spinless" [42l | Fermi particles in N orbitals in the pres- 
ence of fermion-fermion interaction. We shall use second 
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quantization notation with fermion creation and annihi- 
lation operators c\ and fulfilhng the usual anticommu- 
tation relations and creating or annihilating a fermion 
in the ith eigenstate of the single particle Hamiltonian 
defined by the mean field. The M particle eigenstate 
la) = (cj,)"" • • • (c]v_J""-i|0) (where exactly M binary 
digits Qfi of integer a are equal to 1) of the independent 
particle Hamiltonian is then a product of M creation op- 
erators with distinct indices applied to a vacuum state. 
The Hamiltonian is written in the usual manner as 



H = ^ eic\ci + A ^ Vijkic\c]ciCk, 



(4) 



i<j,k<l 



where all Latin indices are running from to — 1, and 
where the are ordered single-particle energies, Vijki 
are properly antisymmetrized two-body matrix elements 
and A is the perturbation strength. The unperturbed 
Hamiltonian Hq will correspond to the Mean Field Ap- 
proximation and will hence be chosen as the one-body 
terms eiC^Q plus the diagonal part of the interaction 

Si<j ^iycJcjcjCi. A natural question is how the re- 
maining part of V , which we shall consider as perturba- 
tion, affects the dynamics. 

The inclusion of the diagonal part of the interaction 
in the unperturbed system is crucial to our argument. 
On one hand it does not affect the eigenfunctions of the 
one-particle term and thus it seems of little importance. 
Yet in the calculation of fidelity the dephasing it pro- 
duces, becomes the dominant term. On the other hand, 
if included in the unperturbed part it will only enter the 
result to the fourth order in A because the spectral two- 
point function will only enter to order A^ and it, in turn, 
will only be affected to order A^ by the diagonal terms. 

So far we have not put any constraints on the spectral 
properties of Hq nor on the perturbation V . The latter 
consists of independent two-body matrix elements, where 
the weight coefficients Vijki are chosen as independent 
random Gaussian complex numbers with vanishing mean 



and {VijkiVr 



ijkl ^ ranopj V 



^ij,opOkl^mn : 



where the variance 



is set to normalize the perturbation as (tr[l/^]) = Af 
where J\f = (^) is the dimension of full Hilbert space. 
The unperturbed dynamics is given by N single-particle 
levels Cj. We shall often choose them as the eigenvalues 
of the NxN Hermitian matrix chosen from the GUE, but 
we shall mention other options and their consequences. 
Similar models have recently been widely studied e.g. in 
nuclear physics [111 and for studying chaotic quantum 
dots [111. 

Two-body operators clcjCiCk can be split into three 
groups by inspecting indices i, j,k,l: diagonal terms with 
two pairs of equal indices, three-orbital terms with one 
pair of equal indices, and four-orbital terms without 
pairs of equal indices. The non-zero matrix elements 
Va/3 in the full Hilbert space therefore couple many- 
particle states which differ by at most two single-particle 
states and these states are the only intermediate states 
in a two- fold transition described by Vj {ti)Vi {t2) in 



After averaging over an ensemble of two-body interac- 
tions all products of independent Vijki vanish and our 
interest is narrowed to transitions with the same initial 
and final state described by averaged matrix elements 
{[V{ti)V{t2)]af:i)v = SapCaih - which only depend 
on the time difference and are denoted by the correlation 
Junction Cair). Using Hq = diag{ha), which is either 



single particle Hamiltonian ho 



cjci, or the mean field Hamiltonian ha ~ (a| '^^hiCi + 
^ Sj<j ^jij^i'^j|Q^)i the correlation function can be writ- 
ten as 

Ca(r) = ^(e2-''"K.^e-2-'«-l/^„)y 



(5) 



The first average is nontrivial in the mean field case 
where the energy differences h-^ — ha contain also di- 
agonal two-body terms (mean field). The correction to 
g-27rii:,((7|ft, l7>-("|ft,|a»e, ^o mean field is of order 
©(A^). As mentioned above, since the correlation func- 
tion will enter the linear response formula for fidelity (26j 
with a prefactor A^, the overall correction to the fidelity 
amplitude is ©(A^) and will be neglected. We split the 
Ca{T) into three groups according to the classification of 
two-body transitions 

Cq(t) = {a\v^^ c\c\ciCkc\c\ciCk 

kl 

+ E 4c]c,C,clc]c,Cfee-2-He.-e.)r (g) 

+ ^ ci.c3'cjC,4c]QCfce-2'^'('='+^^"-^'=-^')na)- 

i<j^k<l 

The abbreviation in the last sum means that indices 
i,j,k,l are all different. The first term of © is actu- 
ally only present if we would be interested in considering 
the full two-body term as perturbation and thus consider 
the corresponding dephasing [i^ . Partially summing the 
expression yields 



Cair) = V '^^^ ] + v 



^(M-l)^nfen,e-2-('=--^'')^ 



+ t--'i^n.n,nfcnie-2-'(^'+^-^'=-^')^ (7) 

ijkl 

Here rn = {a\clci\a) indicates the Fermion occupation 
number and we use the convention = l — rii. The re- 
lation between the fidelity amplitude and the correlation 
function after averaging over the interaction reads 



(/(i)>V' = El^"l'fl-47r2A2 / dt, dt2Caiti-t2) 



(8) 

If the spectrum of Hq is known, this allows us to deter- 
mine the fidelity amplitude by integrating the correlation 
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function in ([6]) twice for each eigenstate a and evaluate 
the sum in ([5]). Actually we are more interested in en- 
semble behavior which we shall study in the next two 
subsections. 



IV. ENSEMBLE AVERAGED FIDELITY 
AMPLITUDE 

Relation ([8]) can be simplified if we average over initial 
states or select a random initial state by choosing Gaus- 
sian random coefficients Xa normalized by the require- 
ment (|a;o,p)o = 1/Af. The correlation function can be 
averaged as C(t) = a (''")) o which impHes aver- 
aging over all distributions of M particles on N orbitals. 
We find 

C{r) = ^(M-l)(f,:') E (°-^"^^"-^"^^)o + 



/iv_4 
~jr[M -2 



(9) 



A closer inspection of first term in ([9]) shows that it is 
connected with the probability distribution of pairs of en- 
ergies of two arbitrary orbitals which is given by Dyson's 
two-point correlation function R2 following the definition 
in [4l| . The average can then be obtained by integration. 
A similar argument connects the second average with 
the four-point correlation function i?4. Finally, using the 
normalization factor = (^/) (2(iV - M) + (^2^^)) we 
rewrite the correlation function in compact form (using 



C{t) = 



+ 




(10) 



The functions !F[Rn] (t) are obtained by Fourier transfor- 
mation of the correlation functions i?„ (see Appendix for 
precise definitions). In the case of GUE single-particle 
spectra, the functions J^[Rn]{t) can be calculated analyt- 
ically and are given as products of finite polynomials in 
t and Gaussians (see Appendix). They are normalized 
such that J^[R2,4] {0) = 1. 

The simulations have been performed both including 
and excluding the diagonal terms of the two-body in- 
teraction and coincide with the theoretical prediction 
in pn]) (or the corresponding full calculation) if the ran- 
dom states are chosen from the whole energy spectrum 
or from its center. The ensemble and state averaged fi- 
delity amplitude can now be obtained by integrating the 
averaged correlation function in (jlOp twice 

{fit)) = 1 - Att^X" f dti / ' dt2C{ti - t2). (11) 
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FIG. 1: Fidelity amplitude decay for 6 particles in 12 orbitals 
averaged over 2 • 10^ realizations of V, Ho and initial states, 
with the perturbation strength A = 10~*. For sufficiently long 
times the decay is quadratic in the case of full perturbation 
(upper) and linear in the residual case (lower). Initial states 
are taken from the center of spectra. 



In the residual case where the perturbation V has no 
diagonal terms (i.e. in the case of mean field unper- 
turbed Hamiltonian) the correlation function has no con- 
stant term. The correlation integral can then be approx- 
imated by a linear function for fairly long times as long 
as the linear response formalism is justified. The sim- 
ulation (Fig. [1]) is done with random initial states from 
the center of the spectra and is in agreement with the 
full-spectra prediction. The residual case thus coincides 
qualitatively with the result of a general residual pertur- 
bation of a Hamiltonian with a random spectrum [iqI . [soj 
also known as POE Fidelity freeze is not found and 
we could easily conclude that this exercise is rather dis- 
appointing. Yet, considering the known non-ergodicity 
of the TBRE, it seems worthwhile to check whether the 
average behavior reflects the typical behavior. 



V. FIDELITY DECAY FOR A TYPICAL 
ENSEMBLE MEMBER 

We therefore proceed to analyze the median of fidelity 
decay. For this purpose we use fidelity F{t) = |/(i)P, 
which is a real quantity, instead of fidelity amplitude 
f{t). (We note however that the results are practically 
the same for the real part of fidelity amplitude.) We de- 
fine median fidelity (t) such that at any time i^m (t) is 
lower than fidelity F{t) in half of the realizations. For a 
randomly chosen member of the ensemble there is thus a 
0.5 probability that a plateau in fidelity decay will last 
longer than the median fidelity plateau (Fig. [2]) . 

A useful alternative in such situation is to consider the 
average of the logarithm of the quantity under consid- 
eration - in our case (ln(l — F{t))). This quantity has 
been used in elastodynamics under the name distortion 
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FIG. 2: Average fidelity in 10^ reafizations of Hq and V 
compared to the median fidelity for 6 particles in 12 orbitals 
and with two values of perturbation strengths A — 10~* and 
A = 10~^. The median fidelity decay (long-dashed, dotted) 
shows a plateau practically absent in the ensemble average 
(full, short-dashed curves). 
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FIG. 3: Fidelity decay in a few individual ensemble mem- 
bers for the residual case with 6 particles and 12 orbitals. 
The upper left inset shows a small region of the correspond- 
ing cumulative level spacing distribution of the unperturbed 
spectra Hq (We plot a relative number of level spacing njM 
below Sn.) The agreement with the minimum spacing only 
theory (dotted) is shown in the upper right inset, dashed line 
denotes the ending time (eq. p6p . 



whereas the lower two curves truly represent the majority 
of cases. The upper left inset of Fig[3] shows the begin- 
ning of the cumulative level spacing distribution of 
corresponding to the fidelity curves shown in the main 
figure. Interestingly, the plateau properties of the excep- 
tional large deviation cases with high plateaus are de- 
termined by the minimum energy spacing alone and the 
conjecture is that the plateau vanishes together with the 
smallest level spacing. 

Another peculiarity results from the two-body nature 
of interaction which only connects certain levels as shown 
in section IIVI Transitions involving three orbitals will 
very likely have big energy difference because of the level 
repulsion in single-particle spectra; this contrasts with 
the transitions involving four orbitals that can easily have 
small energy differences if one particle lowers the energy 
by roughly the same amount by which the other one 
raises it. The effect of a such two-body operator will 
be greater than the effect of any other operator. 

Considering the established connection between the 
level spacing of four-orbital two-body operators and the 
fidelity amplitude plateau, an illustration is made where 
all three-orbital terms are omitted. Each four-orbital 
two-body operator cJc|c/Cfc is in a unique way connected 
with the energy difference [liju :— ei + ej — ek~ei in the 
spectrum of Hq and the four-orbital part of the correla- 
tion function ^ averaged over the perturbation becomes 
a sum over all possible positive spacings 

Ca(t) = v'^ [(a|^znjnznfee-2-l'^-'''l|a) 



(a|nfen;njnie2'''l^'^"'='l|a) 



(12) 



Assuming random initial states averaging can be per- 
formed 



IVc (t) V 



cos2nfiqt (13) 



and the fidelity amplitude is obtained by integration (111!) 



/(t)^l-A^,5^'f ^^l-cos^2.;V^ (14) 



NI{1 



as we have discussed in the introduction. The distortion 
shown in figure 2] in the following section has essentially 
the same feature as the median. Note that it is not ob- 
vious how to obtain analytical results for either of these 
quantities, but the latter may be slightly more accessible. 

It remains to be understood how the freeze, which is 
present in most realizations, comes about and is even- 
tually completely averaged out. To elucidate this prob- 
lem fidelity of a few individual ensemble members is pre- 
sented (Fig. However, the cases shown in the figure 
are not all equally probable; the upper one occurs rarely 



If the smallest degenerate spacing is orders of mag- 
nitude smaller than any other spacing, the sum above 
can be approximated by the smallest spacing term alone. 
All other terms are smaller by a factor (/io/Mg)^- This 
approximation is very illustrative and can be used to es- 
timate the beginning time % and the shift of the plateau 
by the condition /io^b = 1 and by removing the time- 
dependent part, respectively. Estimation of the ending 
time te of the plateau is a little more tedious. There we 
need the fourth-order terms in the linear response for- 
mula (21). Following the same principles and keeping the 
minimum spacing term only we obtain the fourth-order 
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correction to fidelity of highest order in time (we are in- 
terested in very long times) as 



128n^{N -Ay.{N-M-l) XH^ 
N\{N-M)M{M-l){N-M + ?,Y^ 



(15) 



The ending time to can then be estimated by comparing 
amplitudes of the second and fourth-order terms which 
gives the well known [ssl [33| dependence 



^ 2^ v (T' " ^^^^^ " + ^^^^ 



The ending time does not depend on level spacings, which 
is again in agreement with the data shown in Fig. O If 
th ^ tc the plateau, which is a pure second-order phe- 
nomenon, would begin in the region where the second 
order approximation is no longer valid and hence cannot 
exist. The upper right inset on Fig. [3] shows that the 
second and fourth order terms together describe all im- 
portant features of the fidelity amplitude decay for the 
uppermost realization in Fig [3] with an extremely small 
level spacing. 

We now understand why the plateau, nearly always 
present in individual realizations, does not appear in the 
ensemble average. Realizations of the unperturbed spec- 
trum with an extremely small level spacing occur rarely 
but when averaging the fidelity amplitude they eventu- 
ally dominate the behavior and the plateaux are averaged 
out. 



The first model under consideration (uppermost) has a 
random spectrum of Hq such that the level spacing distri- 
bution is Poissonian. The perturbation is chosen from an 
ensemble of full antisymmetric matrices. The Poissonian 
level spacing distribution favors small spacings between 
spectral levels which contribute to fast (initially linear) 
fidelity decay. In this case the type of perturbation has 
no significant effect as small spacings are probable for 
any pair of levels. This is not true for Hq constructed 
from independent single particle spectra which we take 
as Poissonian or GUE as in our original model For 
the latter, small level spacing for some transitions (e.g. 
three-orbital two-body transitions) is very improbable. 
If the perturbation is a full matrix, the effect of such 
transitions is relatively small and there is no significant 
difference between the two cases. If, on the other hand, 
only a few levels are coupled as in the case of two-body 
interactions, three-orbital transitions in the case of Pois- 
sonian single particle spacings can and in the case of 
GUE spacings cannot involve small spacings. Thus in 
the former case they cannot be neglected. Typical decay 
of fidelity is in both cases slower than in the case of full 
perturbation. For illustration, also a completely differ- 
ent physical situation with random matrix many-body 
spectrum of Hq is shown. Because of the level repulsion 
in many-body spectra, any transition will very unlikely 
have small level spacing. Only then, the plateau in fi- 
delity decay is present not only in a typical case but also 
on average. 



VII. GROUND STATE FIDELITY AMPLITUDE 



VI. COMPARISON TO OTHER MODELS 

The theory we have developed to describe the behav- 
ior of the model under consideration can be applied to 
other physical models and used to understand the dif- 
ferences between them. The key features of the model 
are 1) level repulsion in single particle spectra and 
2) sparse perturbation. By comparing to models with 
full random matrix spectra which correspond to chaotic 
dynamics on the one hand, and independent Poissonian 
spectra which corresponds to regular dynamics on the 
other hand, we find that our model lies somewhere in 
between as we can see from the Fig. 2] Besides the me- 
dian fidelity which corresponds to the typical behavior, 
an alternative measure namely the average of the loga- 
rithm of 1 — F{t) (also known a distortion) is considered 
in parallel. Although not completely identical, they both 
show essentially the same features and to some extent 
even quantitative agreement. In order to eliminate any 
potential effect of diagonal perturbation terms, we con- 
sider in the following models always purely imaginary an- 
tisymmetric perturbation, individually normalized in the 
usual way (tr[T/^] = Af). However, we will be switching 
between the cases of a full and a two-body perturbation 
matrix. 



We have up to now considered only random initial 
states, but the behavior of the ground state of the in- 
dependent particle model Hq, i.e. the Hartree-Fock 
ground state, under a perturbation formed by residual 
interactions is also of interest. The correlation func- 
tion can now no longer be simplified by state averaging. 
The main problem is to determine the matrix element 
[Vmmas = Cg{t) where \g) = 4---4_jO) is the 
ground state occupying the lowest M single-particle lev- 
els. Regardless of the realization of the spectra of Hq, 
any non-diagonal two-body operator raises the energy for 
at least the minimum level spacing in the single-particle 
spectra where levels are unlikely to be close together. The 
same incidentally holds in the first excited state with the 
only difference that there the energy can also be lowered 
by the same amount. Note that more freedom does exist 
for higher excited states such as for example 

|s> = 4---4/-24/4f+ilo) (17) 

where the operator c\f_^c\,j_^2'^M+iCM or its hcrmitian 
conjugate increases or decreases energy by eA/_i+eM+2 — 
CM-eM+i = {eM+2-eM+i)-{eM~eM^i) which indeed 
can be very small. 

This leads to the suspicion that the freeze might exist 
on average for these two lowest independent particle ini- 
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FIG. 4: Median fidelity decay (left) and average natural logarithm of 1 — F{t) (right) for various models (described in text 
and labeled in the figure - ordered from top to bottom) for 5 particles in 10 orbitals and perturbation strength A = 10~*. 
Perturbation is always imaginary antisymmetric. The results are averaged over 2 ■ 10^ realizations. 
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FIG. 5: Fidelity amplitude decay for the ground state and 
first excited state compared to higher excited state |s) (eq. 
(|17|) ) for 6 particles in 12 orbitals and perturbation strength 
A — 10^''. The results are averaged over 10'' realizations of 
perturbation and spectra. 

tial states. The numerical results indeed show (Fig. (5]) 
the plateau in the ensemble averaged fidelity amplitude 
for the ground and the first excited state. 



VIII. CONCLUSION 

We have analyzed a very elementary many-fermion 
model in context of echo dynamics. The model de- 
scribes spinless fermions whose underlying single-particle 
dynamics is chaotic and are perturbed by random two- 
body interactions. Our interest was focused on the sta- 
bility of a mean field approximation where the diagonal 
terms of the interaction have been included in the unper- 
turbed Hamiltonian. For weak perturbations the decay 
of the fidelity amplitude in this case typically displays a 



high-level plateau also known as freeze of fidelity. This 
freeze lasts for times long on the scale of the Heisenberg 
time. The unexpected point is that the freeze typically 
present in most realizations of members of the ensemble 
will vanish on average giving a very dramatic example of 
the non-ergodicity of the TBRE. To see the effect which 
should dominate most experiments and should hence be 
observable, we have to consider median behavior or con- 
sider logarithmic averages. This fact beyond the specific 
interest for the TBRE could pave the road for interest- 
ing new applications of random matrix theory analyzing 
non-ergodic situations. 
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APPENDIX: ANALYTICAL CALCULATION OF 
THE CORRELATION FUNCTION 

The final compact expression for the averaged correla- 
tion function C (r) (see eq. [TUl) involves Fourier integrals 
of the Dyson's correlation fmictions i?„(x) 

^[i?„](fc) = I^y^l! |°°i?„(x)e-''=5:?r„H-iH.,d"x. 

(A.l) 

If the correlations functions i?„ correspond to GUE, such 
integrals can be analytically integrated and presented in 
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terms of functions Jnm(t): 

/•oo 

Jnm{t) - / i/„(2;)c-"'i7„(x)e-'*Mx (A.2) 



where Hn{x) are standard Hermite polynomials. It can 
easily be shown using recursion that a function Jrnn{t) 
can be written as finite series 



J mn (^) 



V2"+™n!m! 



n.(m,n) 



E (-2) 



(A.3) 

The simplest case is the transformation of Ri{x) = 

^N-l 1 



Sn=o 2"n!v^ ^ Hn{x)'^ which givcs a simple result 

= Jnn{t). (A.4) 



Higher level correlation functions i?„ are expressible by 
lower level i?„ and cluster functions T2 [4l| . Because T is 
a linear transformation it can also be expressed by lower 
order transformations where transformations of cluster 



functions involve terms without the sign factors (—!)■' in 
the exponential in (|A.1|) . denoted by J^[r„] + . Straight- 
forward calculation gives 

T[T2\^{t) = J2jnm{t)Jnm{Tt) (A.5) 

nm 

nmk 

J'[Ti]^it) = ^ Jml{t)Jrnn{t)Jkl{-t)Jnk{-t). (A.7) 
mnkl 

Finally, following the equality R2{x,y) — Ri{x)Ri{y) — 
T2{x,y) and similarly for i?4 it finally holds 

T[R2]{t) - T[Ri]itf -T[T2]-{t) (A.8) 

T[R,]{t) = T[R,]{t)^+T[T2] + itf+2T[T2]-{tf 

- 2T[Ri]{tfT[T2] + it) - 4T[Ri]{tfT[T2]-{t) 

+ 8T[Ri]{t)T[n]{t)-6T[Ti]it). (A.9) 

Hence, the averaged correlation function, although not 
as a compact expression, can be analytically evaluated 
using the above equations. 
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